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We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse back- 
ground, updating and extending our previous results. A "minimal cosmic reionization model" is 
presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of 
star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases 
rapidly with redshift, from 1% at z = 3 to 50% at z = 8.5, the clumping factor of the intergalactic 
medium is Cigm — 3, and Population III stars make a negligible contribution to the metagalactic 
flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement 
measurements, and yields an optical depth to Thomson scattering, Tg = 0.085, that is in agreement 
with WMAP results. 
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1. Introduction 

The reionization of the all-pervading intergalactic medium (IGM) is a landmark event in the 
history of cosmological structure formation. Studies of Gunn-Peterson absorption in the spectra of 
distant quasars show that hydrogen is highly photoionized out to redshift z > 6, while polarization 
data from the Wilkinson Microwave Anisotropy Probe (WMAP) constrain the redshift of a sudden 
reionization event to be significantly higher, z= ll.Oit 1.4. It is generally thought that the IGM is 
kept ionized by the integrated UV emission from active nuclei and star-forming galaxies, but the 
relative contributions of these sources as a function of epoch are poorly known. In this contribution 
we describe a new version of the code CUBA, aimed at solving the radiative transfer of ionizing 
background radiation in a clumpy, expanding medium. We will use it to compute improved synthe- 
sis models of the UV/X-ray cosmic background spectrum and evolution, combining, updating, and 
extending many of our previous results in this field. For details and relevant references see Haardt 
&Madau (2011). 



2. Cosmological radiative transfer 

The equation of cosmological radiative transfer describing the time evolution of the space- and 
angle-averaged monochromatic intensity Jy is 

d d\ c 

- vH— Jy + 3HJv = -cKyJv + — £v, (2.1) 
ot ov J 471 

where H{z) is the Hubble parameter, c the speed of the light, Ky is the absorption coefficient, and 



Cv the proper volume emissivity. The integration of equation (2.1) gives the background intensity 



at the observed frequency Vo, as seen by an observer at redshift Zo, 



f r \dt/dz\dz^-jP4^ 

AkJz, (l+z)3 



Jyoizo) = jzl \dt/dz\dz\^-^Ey{z)e-\ (2.2) 



where V = Vo(l +z)/(l+Zo), [<i?/<izl =// ^{l+z) ^andT = — ln(e is the effective absorption 
optical depth of a clumpy IGM. 

2.1 IGM absorption 

The effective opacity of the IGM has traditionally been one of the main uncertainties affecting 
calculations of the UV background. Our improved model uses a piecewise power-law parameteri- 
zation for the distribution of absorbers along the hne of sight, 

f{Nm,z)=AN^^{l+z)'', (2.3) 
and is designed to reproduce accurately a number of recent observations. 

• Over the column density range 10" <AfHi< 10'^ cm ^, we use (A, j3, 7) = (1.2 x 10^, 1.5,3.0), 
where the normalization A is expressed in units of cm^^^'^^^^ and j8 = 1.5. 
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• At the other end of the column density distribution, with a power-law exponent J3 = 2 
down to a break column of A'^hi = 10^'^^ cm^, and with an incidence per unit redshift 
°^ (1 the parameters for the dumped Lya systems (DLAs) becomes (A,j3,y) = 
(8.7 X 10l^2, 1.27). 

• For absorbers with 10^^ < A'^hi < 10^*^'^ cm~^ (the so-called "super Lyman-Umit systems", or 
SLLSs), we use (A,j3,7) = (0.45,1.05,1.27). 

• There is obviously a significant mismatch between the power-law exponent for the Lya 

clouds (7= 3) and the SLLSs (7= 1.27). Continuity then requires the shape of f{Nm,z) 
to change with redshift over the colum density range of the Lyman-limit systems (LLSs), 
cm . In this interval of column densities we match the distribution 
function with a power law of redshift-dependent slope. The procedure yields the slopes 
P = 0.47,0.61,0.72,0.82 at redshifts z = 2,3,4,5, respectively. 

• The above parameterizations reproduce well the observations at 2 < z < 5. At low redshift, 
however, HST data show that the forest undergoes a much slower evolution. We take 7=0.16 
in the interval < z < ziow and dN/dz = 34.7 at z = above an equivalent width of 0.24 
A. We derive (A,/3,7) = (1.73 x 10^ 1.5,0.16) for lO" < A^hi < 10^^ cm'^ and (A,j8,7) = 
(5.49 X 10'^2,0.16) for lO^^ < Nui < lO^^-^ cm^^ at all redshifts below ziow = 1.56. We use 
a broken power-law for the redshift distribution of the SLLS s and DLAs as well; assuming 
that the same 7 = 0.16 slope and transition redshift ziow inferred for the forest also hold 
in the case of the thicker absorbers, we derive a normaUzation at z < ziow of A = 1 .28 for 
the SLLSs and A = 2.47 x 10^^ for the DLAs. This yields dN/dz = 0.74 absorbers above 
AThi = 10i^-2 cm-2 at (z) = 0.69. 

• Above z = 5.5 we assume for the forest the values (A, jS, 7) = (29.5, 1.5,9.9) (10^^ <A^hi< 
10^^ cm-2) and (A,j3,7) = (9.35 x 10^2,9.9) (10^^ < A^m < lO^^'^ cm'^) above redshift 
5.5. 

2.2 Source emissivity 

The emissivity is due to several contributing terms: 

• The background photons absorbed through a Lyman series resonance cause a radiative cas- 
cade that ultimately terminates either in a Lya photon or in two-photon 2s ^ Is continuum 
decay. We use the detailed photoionization structure of absorbing systems to calculate the 
reprocessing of background LyC radiation by the clumpy IGM via atomic recombination 
processes. We include recombinations from the continuum to the ground state of HI, Hel, 
and Hell, as well as Hell Bahner, two-photon, and Lya emission. 

• The adopted quasar comoving emissivity at 1 Ryd, £912 (z) /(I + z)^, is 

which closely fits the observational results in the range 1 < z < 5.7 under the assumption 
of pure luminosity evolution. The poorly known faint-end slope of the quasar luminosity 



3 



UV background 



Francesco Haardt 



function at high redshift, incompleteness corrections, as well as the uncertain spectral energy 
distribution (SED) in the UV, all contribute to the large apparent discrepancies between dif- 
ferent authors. We use the functional form given above together with the broken power-law 
SED withLv °< v-O'^'^ for 1300 < A < 5000, andLy °< v'^-^^ for A < 1300. 

• Star-forming galaxies are expected to play a dominant role as sources of hydrogen-ionizing 
radiation at z > 3 as the quasar population declines with lookback time. To compute the LyC 
emissivity from galaxies at all epochs, we start with an empirical determination of the star 
formation history of the universe. We adopt the far-UV (FUV, 1500 A) published luminosity 
functions in the redshift range < z < 9, integrated down to L mm = 0.01 L,, using Schechter 
function fits with parameters {<j>*,L^, cc) to compute the dust-reddened galaxy FUV luminos- 
ity density ppuv- Dust attenuation is treated using a Calzetti extinction law normalized at 
1500 A. Finally, the dust-corrected luminosity densities are smoothed with an approximating 
function and then compared with the results of spectral population synthesis models provided 
by the GALEXEV library 

3. The Spectrum of the UV background 

Figure 1 shows the quasar-only background spectrum generated by an upgraded version of 
our radiative transfer code CUBA, using the formalism and parameters described above. CUBA 
solves the radiative transfer equation by iteration, as its right-hand term implicitly contains J in 
the recombination emissivity and in the effective helium opacity. Physically, this simply means 
that the metagalactic UV flux depends on the ionization state of intervening absorbers, which is in 
turn determined by background radiation. For comparison, we have also plotted the background 
spectrum from our old models. The new models are characterized by a lower UV flux (by as much 
as a factor of 3 at 1 Ryd and z = 3), smaller spectral breaks from HI and Hell LyC absorption, a 
sawtooth modulation by the Lyman series of HI and Hell that becomes more and more substantial 
with redshift. 

Figure 2 shows the full background (i.e., quasar plus galaxies) compared to the quasar-only 
spectrum. The much softer emissivity, linked to HI absorption fixed at the observed level, causes a 
much deeper Hell absorption trough, and much deeper Hell Ly series absorption lines. 

Figure 3 shows our predicted HI ionization rate compared to several measurements. The full 
model provides a good fit to the hydrogen-ionization rates inferred from flux decrement measure- 
ments, predicts that cosmological HE regions overlap at redshift 7.5, and yields an optical depth to 
Thomson scattering, Tg = 0.085, that is in agreement with WMAP results. 
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Figure 1: The cosmic background spectrum from quasars only between 5 A and 5,000 A at epochs z — 
0, 1,3, and 5. The new models (black curves) are compared with the old results of Haardt & Madau (1996, 
red curves). The intensity Jy is expressed in units of 10^^^ ergscm^^s^' Hz^' sr^'. 
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Figure 2: Same as above, but now the red line is the new quasar-only background, compared to the full 
(quasars plus galaxies) model {black line). 
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Figure 3: Hydrogen photoionization rate as a function of redshift. Solid curve: our full model. Dashed 
curves: separate contributions from quasars (leftmost curve) and galaxies (rightmost curve). Different data 
points refer to various empirical measurements from the Lya forest effective opacity (see Haardt & Madau 
20 11 for full references). 
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